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Abstract To develop a targeted metagenomics approach for 
the analysis of novel multispecies microbial consortia in¬ 
volved in the bioconversion of lignocellulose and furanic 
compounds, we applied replicated sequential batch aerobic 
enrichment cultures with either pretreated or untreated wheat 
straw as the sources of carbon and energy. After each transfer, 
exponential growth of bacteria was detected using microscop¬ 
ic cell counts, indicating that the substrate was being utilized. 
In batch, the final bacterial abundances increased from an 
estimated 5 to 8.7-9.5 log 16S rRNA gene copy numbers/ 
ml. The abundances of fungal propagules showed greater 
variation, i.e., between 5.4 and 8.0 log ITS1 copies/ml. Dena¬ 
turing gradient gel electrophoresis analyses showed that the 
bacterial consortia in both treatments reached approximate 
structural stability after six transfers. Moreover, the structures 
of the fungal communities were strongly influenced by sub¬ 
strate treatment. A total of 124 bacterial strains were isolated 
from the two types of enrichment cultures. The most abundant 
strains were affiliated with the genera Raoultella!Klebsiella , 
Kluyvera , Citrobacter , Enterobacter , Pseudomonas , 
Acinetobacter , Flavobacterium and Arthrobacter . Totals of 
43 and 11 strains obtained from the untreated and pretreated 
substrates, respectively, showed (hemi)cellulolytic activity 
(CMC-ase and xylanase), whereas 96 strains were capable of 
growth in 7.5 mM 5-hydroxymethylfurfural. About 50 % of 
the latter showed extracellular oxidoreductase activity as detect¬ 
ed by a novel iodide oxidation method. Also, (hemi)cellulolytic 
fungal strains related to Coniochaeta , Plectosphaerella and 
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Penicillium were isolated. One Tnchospown strain was isolated 
from pretreated wheat straw. The two novel bacterial-fungal 
consortia are starting points for lignocellulose degradation 
applications. 

Keywords Enrichment cultures • Enterobacteriales • 
Lignocellulose • 5-HMF • Oxidoreductase activity 

Introduction 

Industries around the world have been recently focusing on 
improvements in the production of, for instance, fuel (etha¬ 
nol-oil-gas) from biological waste materials like lignocellu- 
losic compounds (Xing et al. 2012). However, processes are 
often suboptimal and several key issues need attention, like 
the high cost of thermochemical conversion, the efficiency of 
enzymatic conversion and the toxicity of some byproducts 
that emerge in pretreatments of source material (Zheng et al. 
2009). Lignocellulosic plant waste materials are composed of 
cellulose, hemicellulose, lignin and pectin, constituting rather 
inaccessible structures. Thus, substrate pretreatment is thought 
to be required to open up and/or remove the lignin moieties, 
depolymerize cellulose and hemicellulose and make the ma¬ 
terial amenable to enzymatic attack. Unfortunately, pretreat¬ 
ment steps often generate byproducts such as furfural, 5- 
hydroxymethylfurfural (5-HMF) and phenolic compounds 
(vanillin, syringaldehyde and coniferyl aldehyde), which are 
inhibitors of subsequent transformation steps (Chandel and 
Singh 2011). 

The microbial conversion of lignocellulosic substrates to 
monomers involves the synergistic action of several enzymes, 
such as peroxidases, laccases, xylanases, endoglucanases, 
exoglucanases and (3-glucosidases (Hasunuma et al. 2013). 
In addition, furfural/5-HMF oxidoreductase, 2,5-furan dicar- 
boxylic acid decarboxylase, 2-furoyl-CoA synthethase, 
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fiiroyl-CoA dehydrogenase and 2-oxoglutaryl-CoA hydrolase 
are involved in conversions of furfural and 5-HMF to 2- 
oxoglutaric acid, which is metabolized via the tricarboxylic 
acid pathway (Wierckx et al. 2011). In recent years, a high 
demand for efficient enzymes that allow the bioconversion of 
different lignocellulosic substrates and furanic compounds has 
emerged. Hence, novel genes/enzymes and microorganisms 
capable of metabolizing both types of compounds and 
fermenting sugars are deemed to be very important for use 
in consolidated bioprocessing (CBP) technologies for biofuel 
production (Cheng and Zhu 2012). 

There are different strategies in the search for efficient 
microorganisms that degrade lignocellulosic matter: (a) devel¬ 
opment of specialized microbial communities using enrich¬ 
ment cultures coupled or not to stable isotope probing (SIP) 
and (b) screening for enzymes via functional metagenomics. 
Recently, a concept called habitat biasing was introduced, in 
which specific substrates (in situ or ex situ) would incite the 
emergence of microbial communities fine-tuned to the desired 
function (Ekkers et al. 2012). This methodology has great 
potential, in the sense that microbial consortia with intricate 
links between them, providing stability, are obtained. Yet, a 
possible drawback is that particular members of the commu¬ 
nity that become prevalent are mere scavengers of unlocked 
carbon, i.e., behave like “cheaters” in the process. Different 
approaches to isolation and enrichment of functional microbial 
consortia have recently been reported (Lee et al. 2013). For 
example, some systems in which hemicellulose degraders were 
enriched have been set up with sugarcane bagasse, switchgrass, 
or wheat straw as carbon sources (Wongwilaiwalin et al. 2010; 
Gladden et al. 2011; Wang et al. 2011a). In addition, furanic 
aldehyde degrading microorganisms have been isolated using 
classical enrichment on furfural and/or 5-HMF (Lopez et al. 
2004) and also in torrefied grass biomass (Trifonova et al. 
2008). 

A large number of mesophilic aerobic bacteria are known 
to be capable of degrading the lignocellulosic material present 
in plant cell walls. The degrading organisms generally belong 
to the orders Actinomycetales , Rhizobiales and Bacillales. 
However, recently, the importance of members of the 
Enterobacteriales (i.e., Klebsiella, Enterobacter, Citrobacter 
and Pantoea), Pseudomonadales and Sphingobacteriales has 
been indicated (Rattanachomsri et al. 2011; Aylward et al. 
2012; Wang et al. 2013; Kanokratana et al. 2013). In addition, 
members of the Ascomycota and Basidiomycota have also 
been described as key lignocellulose degraders in soils 
(Stursova et al. 2012). Also, Trichoderma , Aspergillus , Rhi- 
zopus,Mucor , Penicillium and Trametes types have been the 
focus of much previous research. The most abundant furanic 
compound degrading bacteria are members of the orders 
Burkholderiales , Pseudomonadales , Actinomycetales and 
Xanthomonadales (Wierckx et al. 2011), whereas three 
fungal furanic aldehyde degraders have been reported, i.e., 


Coniochaeta ligniaria (Lopez et al. 2004), Trichosporon 
cutaneum (Chen et al. 2009) and Amorphotheca resinae 
(Zhang et al. 2010). 

Given the likely importance of the behavior of the micro¬ 
bial consortia involved, we here posit that in-depth studies on 
these will provide an important foundation for understanding 
the complex interactions of lignocellulosic degradation in 
nature and between bacterial and fungal players, providing a 
platform for biotechnological applications and CBR This 
study aimed to produce novel microbial consortia involved 
in the bioconversion of lignocellulosic substrates and furanic 
toxic compounds. The focus was on the production of stable 
bacterial-fungal consortia that are growing on either 
pretreated or untreated wheat straw substrate. We here de¬ 
scribe their characteristics and the isolation of some of the 
key members. The novel consortia may constitute starting 
points for 16S rRNA gene-ITS pyrosequencing and 
metagenomic analyses. 

Methods 

Substrate preparation 

Wheat straw material was collected from a farm near Gro¬ 
ningen (The Netherlands) and air-dried before cutting into 
pieces of about 5 cm length. Two different substrates were 
thus obtained, to serve as carbon and energy sources: (1) 
“raw” wheat straw (RWS1) and (2) torrefied “mild pyrolysis” 
wheat straw (TWS1). The torrefaction process (using 20 g of 
dry wheat straw) was performed in a carbolite muffle furnace 
at a temperature of 240 °C during 1 h, in accordance with 
Trifonova et al. (2008). Lignocellulosic materials (RWS1 and 
TWS1) were ground, using a hammer mill, to pieces < 1 mm to 
augment the surface capacity access and were added in the 
enrichment medium described below. 

Construction of enriched cultures by dilution-to-stimulation 
approach 

Ten randomly taken soil samples of 10 g were collected and 
mixed from a forest soil (0 to 10 cm depth) in Groningen, The 
Netherlands (53.41 N; 6.90 E) in August 2012 (summer 
season). At the sampling site, a high quantity of decaying 
wood was present (Supplementary Fig. SI). Cell suspensions 
were prepared by adding 10 g of sampled soil to 250 ml flasks 
containing 10 g of sterile gravel in 90 ml of mineral salt 
medium (MSM) (7 g/1 Na 2 HP0 4 ; 2 g/1 K 2 HP0 4 ; 1 g/1 
(NH 4 )S0 4 ; 0.1 g/1 Ca(N0 3 ) 2 ; 0.2 g/1 MgCl 2 , pH 7.5). The 
flasks were shaken for 20 min at 250 rpm and 3 ml of 
suspension was then sampled for DNA extraction. Aliquots 
(250 ql) of soil suspension were added to triplicate Erlenmey- 
er flasks containing 25 ml of MSM with 1 % lignocellulose 
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substrate and supplemented with 25 pi of trace element (2.5 g/1 
EDTA; 1.5 g/1 FeS0 4 ; 0.025 g/1 CoCl 2 ; 0.025 g/1 ZnS0 4 ; 
0.015 g/1 MnCl 2 ; 0.015 g/1 NaMo0 4 ; 0.01 g/1 NiCl 2 ; 0.02 g/1 
H 3 B0 3 ; 0.005 g/1 CuCl 2 ) and 25 jllI of vitamin solution (0.1 g/1 
Ca-panthothenate; 0.1 g/1 cyanocobalamine; 0.1 g/1 nicotinic 
acid; 0.1 g/1 pyridoxal; 0.1 g/1 riboflavin; 0.1 g/1 thiamin; 
0.01 g/1 biotin; 0.1 g/1 folic acid). Flasks were incubated at 
25 °C in aerobic conditions, with shaking at 100 rpm. Two 
controls, i.e., one without substrate (WS) and the other 
without microbial source (WMS), were also set up (Sup¬ 
plementary Fig. SI). Cultures were monitored at regular 
time intervals and once systems reached high cell density 
(7-8 log cells/ml) and straw was visually degraded (be¬ 
tween 6 and 8 days), aliquots (25 \i\ of microbial suspen¬ 
sion with fibrous material) were transferred to 25 ml of 
fresh wheat straw containing medium (dilution 10 -3 ). This 
procedure was performed ten times. The numbers of microbial 
cells per milliliter directly after the transfer and at the end of 
growth in each transfer were quantified by microscopical cell 
counting. Finally, samples were taken from the consortia at 
each transfer and stored with 20 % of glycerol at -80 °C. All 
chemicals and reagents used in this work were of analytical or 
molecular biology grade. 

DNA extraction and real-time PCR quantification 

DNA extractions from the cultures were performed using 3 ml 
samples. The Power Soil DNA extraction kit (MoBio® Fab- 
oratories Inc., Carlsbad, USA) was used according to the 
manufacturer’s instructions. The only modification made was 
the resuspension of extracted DNA in 75 pi The 16S rRNA 
gene (region V5-V6) and ITS1 region were amplified by 
qPCR using 1 ng of diluted extracted DNA as a template 
and primers 16SFP/16SRP and 5.8S/ITS1, respectively 
(Pereira e Silva et al. 2012). Standard curves were constructed 
using serial dilutions of plasmid (1 to 8 log copies/pi) that 
contained cloned 16S rRNA and ITS1 fragments from 
Serratia plymuthica (KF495530) and Coniochaeta ligniaria 
(KF285995), respectively. Absolute quantification was car¬ 
ried out in three replicates on ABI Prism 7300 Cycler (Ap¬ 
plied Biosystem, Fohne, Germany). 

PCR-DGGE analysis 

Denaturing gradient gel electrophoresis (DGGE) was 
performed in the Ingeny Phor-U System (Ingeny International, 
Goes, The Netherlands) following PCR with primers F968 
with a GC clamp attached to the 5' end and the universal 
bacterial primer R1401.1b for 16S rRNA gene. For fungal 
communities, primers EF4/ITS4 were used, followed by a 
second amplification with the primers ITSlf-GC/ITS2. Primer 
sequences, PCR mixtures and cycling conditions were de¬ 
scribed by Pereira e Silva et al. (2012). The DGGE was 


performed in 6 % (w/v) polyacrylamide gels with 45-65 and 
20-50 % denaturant gradient for bacterial and fungal commu¬ 
nities, respectively (100 % is defined as 7.0 M urea with 40 % 
deionized formamide). Electrophoresis was carried out at 
100 V for 16 h at 60 °C and the gels were stained for 
30 min in 0.5 % TAE buffer with SYBR gold (final concen¬ 
tration of 0.5 qg/1) (Invitrogen, Breda, The Netherlands). 
Images were taken using Imagemaster VDS (Amersham Bio¬ 
sciences, Buckinghamshire, UK). Fingerprinting results were 
analyzed using GelCompar software (Applied Maths, Sint- 
Martens Fatem, Belgium) and then principal components 
analysis (PCA) was performed against the DGGE banding 
pattern using Canoco (v4.52; Wageningen, The Netherlands) 
(Andreote et al. 2009). 

Recovery and identification of bacterial and fungal strains 

The isolation of bacterial and fungal strains along with trans¬ 
fers was done on R2A agar (BD Difco®, Detroit, USA) and 
Potato Dextrose agar (PDA) (Duchefa Biochemie BV, Haar¬ 
lem, The Netherlands), respectively. Serial dilutions were 
done in MSM and 100 pi of dilution (10 -7 to 10 -9 ) was spread 
on the surface of each medium. Morphological different col¬ 
onies (Supplementary Fig. SI) were purified and preserved at 
4 and -80 °C (in R2A broth with 20 % of glycerol). The 
genomic DNA was obtained using Microbial DNA Isolation 
Kit (MoBio®). For bacterial strains, we performed a BOX- 
PCR in order to cluster the strains, with the primer AIR 
(CTACGGCAAGGCGACGCTGA). BOX-PCR amplifica¬ 
tion was performed under the following conditions: 95 °C 
for 2 min; 35 cycles of 94 °C for 3 s, 50 °C for 1 min, 65 °C for 
8 min; and a final extension step at 65 °C for 16 min. BOX 
clustering analyses were performed using GelCompar soft¬ 
ware. Bacterial 16S rRNA genes of representative strains that 
belonged to each BOX group were amplified using primers 
27F (5'-AGAGTTTGATCMTGGCTCAG-3') and 1492R 
(5'-GGTTACCTTGTTACGACTT-3'). PCR reactions were 
done in a 50-pi reaction mixture containing 1.25 U of Taq 
DNA polymerase (Biolase®, Fuckenwalde, Germany), lx 
reaction buffer, 2 mM of MgCl 2 , 200 qM of dNTPs, 0.3 \iM 
of each primer, 0.25 pi of formamide (100 %), 0.125 pi of 
BSA (20 mg/ml) and 50-100 ng of DNA. PCR reactions were 
run for 35 cycles as follows: denaturation at 95 °C for 45 s, 
annealing at 54 °C for 1 min and elongation at 72 °C for 2 min. 
For fungal strains, we amplified the partial 18S rRNA region 
using primers EF4 and ITS4 according to Pereira e Silva et al. 
(2012). PCR products were purified with a Wizard SV Gel 
and PCR Clean Up system (Promega, Madison, USA) and 
sequenced by Sanger technology (FGC Genomics, Germa¬ 
ny) using the 27F (for bacteria) and ITS4 (for fungi) 
primers. All chromatograms were analyzed for quality using 
the CFC Main Workbench software (CFCbio, MD, USA). 
Taxonomic assignment of the sequences was done by using 
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BLAST-N against the NCBI Database (http://blast.ncbi.nlm. 
nih.gov/Blast.cgi). 

Screening for (hemi)cellulolytic activity, growth in furanic 
compounds and 5-HMF oxidoreductase activity 

Briefly, a small amount of bacterial pure colony (on R2A agar) 
was transferred to 200 pi of R2A broth using a sterile toothpick. 
Cultures in microplates were grown overnight and shaken at 
100 rpm at 25 °C. Screenings for (hemi)cellulolytic activity and 
growth on fiiranic compounds were done in mineral minimum 
agar (MMA) (0.2 % NaN0 3 ; 0.1 % K 2 HP0 4 ; 0.05 % MgS0 4 ; 
0.05 % KC1; 1 % of vitamin solution; 1.5 % agar). We evaluated 
the growth (negative, weak and positive) of the strains in the 
presence of 0.2 % glucose, 0.2 % fmctose, 0.2 % carboxymeth- 
ylcellulose (CMC), 0.2 % xylan from beechwood, 7.5 mM 
furfural and 7.5 mM 5-HMF as a final concentration. A drop 
(20 pi) of bacterial culture grown overnight was inoculated on 
an agar plate (16 strains/per plate). Fungal strains were grown 
on PDA, after which agar plugs (5 mm in diameter) containing 
grown mycelium were excised from the plate of each isolate. 
The agar plugs were placed in the center of each agar plate 
surface. All assays were performed in duplicate using as a 
negative control MMA without carbon source. The plates were 
incubated at 28 °C for 36 h and, after evaluation of growth, 
were flooded with Gram iodine (Kasana et al. 2008) for detec¬ 
tion of CMC-ase and xylanase activity and evaluation of oxi¬ 
doreductase activity in the presence of 5-HMF. We screened 
totals of 124 and 7 bacterial and fungal strains, respectively. 
CMC and xylan degradation was indicated by detection of 
clearing zones around the colonies and an enzymatic activity 
index (EAI) was defined on the basis of the diameter 
(millimeter) of the zone of clearance plus the diameter of the 
colony divided by the diameter of the colony. A cut-off value of 
more than 2.0 mm was considered as a positive. The same 
index was applied to describe the oxidoreductase activity. 

Nucleotide sequence accession numbers 

Sequences are publicly available in the GenBank database 
under accession numbers KF285950-KF285998. 


Results 

Enrichment cultures and real-time PCR cell quantification 

Microcosm systems were constmcted with cells from a forest 
soil as the inoculum (6 log cells/ml in the first enrichment 
flasks). Bacterial growth was observed after each transfer to 
fresh substrate-containing medium the cell densities increas¬ 
ing from the inoculum level, around 5 to 8.0-8.4 log bacterial 
cells/ml in all transfers and for both strategies (Fig. 1). In the 


^2 TWSl (starting inoculum) TWS1 (achieved growth) I I RWS1 (starting inoculum) I RWS1 (achieved growth) 



Fig. 1 Bacterial cell numbers (log cells per milliliter) in enrichment 
cultures with RWS1 and TWSl. Starting inocula revealed growth after 
6-8 days in every sequential transfer step 


without-substrate negative control (WS) the number of cells 
per milliliter decreased from 5 to 2 log cells/ml from the third 
to the sixth transfers, respectively (data not shown). 

The numbers of 16S rRNA gene and ITS 1 copies in the forest 
soil suspension used as the inoculum were around 8 log per ml 
(1x10 s copies/ml, ending up as lxlO 6 copies/ml in the first 
enrichment flasks). In the first flasks, the 16S rRNA gene copy 
numbers increased to 8.9 (±0.83) and 9.2 (±0.83) log per ml for 
the RWS1 and TWSl, respectively (Fig. 2a), whereas these 
remained reduced (around 7 log copies/ml) in the WS control 
as compared with the enrichment strategies. The abundances of 
bacterial communities consistently ended up at gene copy num¬ 
bers between 8.7 and 9.5 log copies per ml following each 
transfer and subsequent growth (Fig. 2a). An assessment of the 
16S rRNA gene copy numbers per cell in transfer 10 for RSW1 
(4.53 x 10 9 copies/ml and 1.55 x 10 9 cells/ml) revealed an average 
of 2.92 copies of 16S rRNA genes per cell. However, in the last 
transfer for the TWS1 system, this value was 1.71 copies per cell. 

In the first transfer flasks, the fungal communities showed 
numbers of 5.6 (±0.66) and 5.9 (±0.87) log ITS1 copies/ml in 
RWS1 and TWSl, respectively (Fig. 2b), whereas around 6 log 
copies/ml had been added. Fungal communities showed a 
greater variation in abundance, with maxima in transfer-7 (8.0 
±0.26 log copies/ml) and transfer-4 flasks (7.8±0.42 log 
copies/ml) for TWSl and RWS1, respectively. However, low 
values were observed in the flasks of transfer 6 (5.8± 1.32 log 
copies/ml) and 4 (5.4±0.30 log copies/ml) for TWSl and 
RWS1, respectively (Fig. 2b). Moreover, the qPCR assess¬ 
ments showed progressively lower numbers in the control 
flasks (WS, WMS) along with the transfers, indicating that 
low or no growth took place without added wheat substrate or 
microbial source. Thus, every 1:1,000 transfer did effectively 
decrease the number of target gene copies, used as proxies for 
cell numbers, in a progressive manner. 

PCR-DGGE analyses 

We used PCR-DGGE to evaluate bacterial and fungal com¬ 
munity structures, as well as their stabilities, along with all 
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H Soil suspension DrWSI ■ TWS1 □ WMS (RWS1) E WS 



Number of transfers 


H Soil suspension QrWSI ■ TWS1 □ WMS (RWS1) 0 WS 
12 - 


1( H 



Soil a 1 23 45678 10 


Number of transfers 


Fig. 2 Abundance of bacterial and fungal communities as measured by 
qPCR in forest soil—derived enrichment cultures with RWS1, TWS1 and 
microbial source or substrate-less negative controls (respectively, WMS 
and WS). a 16S rRNA gene and b ITS 1 copy numbers (log copy number 
per milliliter of culture) in the sequential transfer steps. Quantifications 


were performed after completion of growth in each transfer, which 
commonly took 6-8 days. Superscripted letter a indicates the quantifi¬ 
cation in the initial soil suspension. Bars refer to standard errors of the 
mean (n =3) 


sequential transfers (Fig. 3). For RWS1, the bacterial commu¬ 
nities revealed low diversity structures as from the beginning, 
which reached stability after six transfers. In the first transfer 
for this system, we observed approximately 12 bands and the 
profiles fluctuated a bit along six transfers, revealing minor 
changes in the number of bands (between 10 and 14). Three 
bands were unique and dominant in the RWS1 profiles after 
the third transfer (Bl, B3 and B5) (Fig. 3a). Then, the profiles 
showed approximately 12 to 14 major bands between transfer 
6 and transfer 10. For TWS1, eight bands dominating the 
profiles were observed in the first, second and third transfers. 
Approximately ten bands were observed after six to ten trans¬ 
fers, with around 50 % of similarity between the third and 
tenth transfers. In the TWS1, one dominant band (B2), possi¬ 
bly belonging to Bacteroidetes , appeared in the second trans¬ 
fer flask and remained dominant throughout the further flasks. 
Five dominant bands (B1-B5) from the two enrichment strat¬ 
egies were retrieved and sequenced. The dominant and 


consistent DGGE band B4 from TWS1 (transfer 6) was 
assigned to Raoultella terrigena (99 % identity— 
NR037085.1); it was also dominant in RWS1 after transfer 
3. Unfortunately, sequencing of the other bands failed. The 
triplicates of each system, in each transfer, consistently re¬ 
vealed communities that were very similar to each other (flask 
to flask) in both strategies (Fig. 3 a). 

With respect to the fungal community structures, the PCR- 
DGGE results showed (1) a high fungal community richness 
in the original soil sample (approximately 30 bands) and (2) 
consistency between replicates (flask to flask) in each transfer 
with reduced richness. In the first transfer for RWS1, multiple 
bands were observed, i.e., 20 to 25 bands, whereas these 
numbers dropped progressively upon transfer. At transfer 8, 
we observed 8 to 12 different bands. For TWS1, between 
eight and ten bands were observed in the first stages, 
suggesting that a few specific fungi had been selected 
(Fig. 3b). The profiles and band numbers were approximately 
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stable along transfer 3 to transfer 8, with six to eight bands in 
the last transfers. Specifically, two large clusters could be 
discerned in the TWS 1-derived profiles from all flasks, which 
shared around 40 % similarity. The first cluster encompassed 
all profiles derived from transfers 1 to 5, whereas the second 
one encompassed all samples from transfers 6 to 8 (Fig. 3b). 


Recovery and identification of bacterial and fungal strains 

Totals of 60 bacterial strains were recovered from the systems of 
transfers 1, 2, 3, 6 and 8 (30 strains for both cultures). From the 
transfer-10 flasks, we recovered 64 bacterial strains (32 strains 
for both cultures). Pre-grouping of all strains by BOX-PCR 
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revealed a total of 13 large groups (data not shown). We then 
selected representatives from each BOX-PCR group and se¬ 
quenced their 16S rRNA genes for identity (Table 1). The 
RWS1 enrichment systems showed high diversities and strains 
were mainly affiliated to species of the genera Raoultella/Kleb¬ 
siella , Kluyvera, Citrobacter, Enterobacter, Serratia 
( Enterobacteriales ), Pseudomonas, Acinetobacter, 
Flavobacterium, Aeromonas and Arthrobacter. From the 
TWS1 system, we identified seven bacterial types affiliated with 
species of Raoultella/Klebsiella, Enterobacter, Citrobacter, 
Escherichia, Pseudomonas, Sphingobacterium and 
Arthrobacter (Fig. 4 and Table 1). The organisms identified as 
BOX-PCR groups 4,6,8,10 and 12 (i Citrobacter, Enterobacter, 
Kluyvera, Raoultella/Klebsiella and Escherichia) also showed 
high similarities with other genera within the order 
Enterobacteriales. However, we used the “best” BLAST hit 
affiliation for taxonomic assignment (Table 1) using a cutoff of 
90 % of identity and 95 % of coverage. It is known that species 
of different genera from the Enterobacteriales are genetically 
diverse, despite the near identity (97-99 %) of the 16S rRNA 
gene sequence and so the classifications are preliminary and will 
need confirmation with data from, e.g., comparative genomic 
hybridization and phenotype analyses. We isolated fungal strains 
from the systems of transfers 2, 8 and 10. In the transfer-2 
system, in both wheat substrates, we obtained strains that were 
closely related to Coniochaeta, Plectosphaerella and 
Trichosporon species, whereas in the transfer-8 and transfer-10 
systems, we identified strains belonging to species of Penicilli- 
um (Fig. 4 and Table 2). 

(Hemi)cellulolytic activities of bacterial and fungal strains 

A total of 124 bacterial and 7 fungal strains was evaluated for 
(hemi)cellulolytic activity (CMC-ase and xylanase). The data 
obtained showed that all fungal strains possessed 
(hemi)cellulolytic activity, with the exception of one isolate 
of Trichosporon (Table 2 and Supplementary Fig. S2). All 
bacterial strains were able to grow (strong to weak) in 0.2 % 
fructose and glucose (Fig. 5a). Totals of 21 (75 %), 6 (75 %), 4 
(80 %) and 4 (100 %) RWS1-derived strains of Raoultella/ 
Klebsiella, Kluyvera, Enterobacter and Flavobacterium , re¬ 
spectively, showed positive CMC-ase and xylanase activity. 
All strain groups retrieved from RWS1 had members with 
(hemi)cellulolytic activity, except for Pseudomonas . Only six 
(19.3 %) and five (29.4 %) of the TWS1-isolated strains 
related to Raoultella/Klebsiella and Pseudomonas , respec¬ 
tively, showed enzymatic activity (Fig. 5b). 

Growth in furanic compounds and oxidoreductase activity 
in the presence of 5-HMF 

All bacterial strains were evaluated for their ability to grow in 
7.5 mM of 5-HMF or furfural agar. After 2 days at 27 °C, some 


showed residual growth on furfural agar similar to that of the 
negative control. Thus, we considered this result as a negative. A 
total of 96 strains showed colonies with white concentric circles 
of growth in 5-HMF (Fig. 5a). From RWS1, 27 (96.4 %), 5 
(100 %) and 5 (71.4 %) strains of Raoultella /Klebsiella, 
Enterobacter and Acinetobacter , respectively, were capable of 
growth in 5-HMF. Only one strain, affiliated with Aeromonas, 
was not able to grow. From TWS1, totals of 28 (90.3 %), 9 
(52.9 %) and 5 (100 %) strains of Raoultella/Klebsiella, Pseu¬ 
domonas and Arthrobacter showed growth in 5-HMF (Fig. 5b). 

Thus, after we observed growth of bacterial strains, we 
tested the reaction of the colonies by Gram stain iodine 
(33 % of I 2 and 67 % of KI) on 5-HMF plates. After flooding 
the colonies, we observed the production of a yellow halo in 
some strains (gray scale color plates in Fig. 5a). Two concepts 
may be put forward to explain this response: (1) the produc¬ 
tion of H 2 0 2 -oxidation of 2KI in acidic conditions-produc- 
tion of I 2 and change of color (Schmitz 2010; Milenkovic and 
Stanisavljev 2011) or (2) production of a 5-HMF oxidoreduc¬ 
tase (HmfH)-direct oxidation of 2KI-production of I 2 and 
change of color (Amachi et al. 2005; Suzuki et al. 2012). In 
both concepts, we expect to detect extracellular oxidoreductase 
activity (Fig. 5). Thus, 32 and 15 bacterial strains, obtained 
from RWS1 and TWS1, respectively, showed extracellular 5- 
HMF oxidoreductase activity (Fig. 5b). Remarkably, we did not 
observe yellow halo formation on MMA with glucose, fructose 
or furfural as a carbon source (data not shown). 

Discussion 

The biodegradation of lignocellulosic biomass and subsequent 
production of green chemistry/renewable fuels through complex 
communities that work together in the presence of toxic com¬ 
pounds has been proposed as a very promising biotechnological 
approach, which, however, is in continuous need of improve¬ 
ment (Zuroff and Curtis 2012). In this study, two different ex situ 
enrichment strategies to produce relevant microbial consortia 
were used with (1) untreated wheat straw (RWS1) and (2) 
torrefied wheat straw (TWS1). Unlike dilution-to-extinction 
approaches, which incur the loss of a certain function from a 
progressively diluted inoculum, the approach used by us (dilu- 
tion-to-stimulation) results—due to selective processes—in a 
stimulus of function within the emerging consortia during suc¬ 
cession (Lee et al. 2013; Zhang et al. 2013). In a previous work, 
microbial communities have been enriched from soil on NaOH- 
treated wheat straw (Wang et al. 2011a), but only a few studies 
have used untreated substrate (Hui et al. 2013). Wheat straw is 
mainly composed of cellulose (35-50 %), hemicellulose (25- 
30 %) and lignin (25-30 %) (Cheng and Zhu 2012), inciting us 
to use torrefaction to loosen the lignin moieties (Nikolopoulos 
et al. 2013). However, this possibly released furanic compounds. 
In a previous work (Trifonova et al. 2008), we applied 
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Table 1 

Closest matches for 16S rRNA gene sequences of strains in each BOX group and relative abundance in RWS1 and TWS1 


Strain 

code 

Box 

group 

Closest relative species—“best” 
hit partial 16S rDNA (accession) 

Length 

(bp) 

Coverage 

(%) 

Identity 

(%) 

Closest relative 
genus 

NSRWS1 
(% relative 
abundance) 

NSTWS1 
(% relative 
abundance) 

10w4 

1 

Acinetobacter baumannii (NR 074737.1) 

964 

99 

96 

Acinetobacter 

7(11.29) 

0 (0.0) 

10wl6 


Acinetobacter calcoaceticus (NR 0423 87.1) 

966 

99 

97 




lw4 


Acinetobacter calcoaceticus (NR 0423 87.1) 

930 

100 

99 




6wl 


Acinetobacter haemolyticus (NR 026207.1) 

902 

99 

95 




8w3 


Acinetobacter johnsonii (NR 044975.1) 

791 

100 

98 




8w5 


Acinetobacter johnsonii (NR 044975.1) 

874 

100 

98 




6w3 

2 

Aeromonas vemnii (NR 044845.1) 

1,095 

100 

99 

Aeromonas 

1 (1.61) 

0 (0.0) 

2t3 

3 

Arthrobacter aurescens (NR 074272.1) 

1,001 

99 

97 

Arthrobacter a 

2 (3.22) 

5 (8.06) 

lw3 

4 

Citrobacter gillenii (NR 041697.1) 

1,044 

100 

99 

Citmbacter 

6 (9.67) 

2 (3.22) 

10t3 

6 

Entembacter amnigenus (NR 024642.1) 

974 

99 

98 

Entembacter 

5 (8.06) 

3 (4.83) 

10w26 


Entembacter amnigenus (NR 024642.1) 

983 

99 

98 




lw6 


Entembacter amnigenus (NR 024642.1) 

883 

100 

98 




lt2 


Entembacter asburiae (NR 074722.1) 

771 

100 

98 




3w2 

7 

Flavobacterium glaciei (NR 043891.1) 

572 

100 

92 

Flavobacterium 

4 (6.45) 

0 (0.0) 

10w3 


Flavobacterium johnsoniae (NR 074455.1) 

973 

99 

93 




6w6 

8 

Kluyvera cryocrescens (NR 028803.1) 

444 

98 

93 

Kluyvera 

8 (12.90) 

0 (0.0) 

10wl8 


Kluyvera intermedia (NR 028802.1) 

987 

100 

98 




10wl9 


Kluyvera intermedia (NR 028802.1) 

988 

99 

98 




10w21 


Kluyvera intermedia (NR 028802.1) 

1,014 

99 

97 




2w3 


Kluyvera intermedia (NR 028802.1) 

916 

100 

98 




6w2 


Kluyvera intermedia (NR 028802.1) 

939 

100 

98 




3t3 

9 

Pseudomonas asplenii (NR 040802.1) 

827 

100 

97 

Pseudomonas 

1 (1.61) 

17 (27.41) 

8t4 


Pseudomonas oryzihabitans (NR 025881.1) 

1,032 

99 

96 




10tl5 


Pseudomonasplecoglossicida (NR 024662.1) 

1,008 

99 

97 




10w8 


Pseudomonas putida (NR 043424.1) 

1,071 

100 

98 




2t2 


Pseudomonas trivialis (NR 028987.1) 

984 

99 

94 




10t20 


Pseudomonas vranovensis (NR 043313.1) 

1,030 

96 

95 




10tl2 

10 

Raoultella terrigena (NR 037085.1) 

996 

100 

99 

Raoultella (Klebsiella-Yike) 

28 (45.16) 

31 (50.0) 

10tl4 


Raoultella terrigena (NR 037085.1) 

1,044 

100 

98 




10tl7 


Raoultella terrigena (NR 037085.1) 

1,010 

99 

98 




10t5 


Raoultella terrigena (NR 037085.1) 

822 

99 

99 




lOwll 


Raoultella terrigena (NR 037085.1) 

990 

99 

98 




2w5 


Raoultella terrigena (NR 037085.1) 

992 

99 

97 




6tl 


Raoultella terrigena (NR 037085.1) 

995 

100 

99 




8wl 


Raoultella terrigena (NR 037085.1) 

1,002 

99 

99 




8w6 


Raoultella terrigena (NR 037085.1) 

1,053 

99 

98 




2w4 

11 

Serratia pmteamaculans (NR 025341.1) 

806 

99 

99 

Serratia 

1 (1.61) 

0 (0.0) 

lt5 

12 

Escherichia fergusonii (NR 074902.1) 

752 

100 

97 

Escherichia 

0 (0.0) 

2 (3.22) 

ltl 


Escherichia fergusonii (NR 074902.1) 

1,065 

99 

98 




3t5 

13 

Sphingobacterium faecium (NR 025537.1) 

856 

99 

97 

Sph ingobacterium 

0 (0.0) 

2 (3.22) 

3t6 


Sphingobacterium faecium (NR 025537.1) 

906 

100 

98 





Superscripted letter a belongs to Actinobacteria 

NSRWS1 number of strains recovered from RWS1, NSTWS1 number of strains recovered from TWS1 


torrefaction to grass fibers and found furfural (27.9 mg/kg) and 
5-HMF (<10 mg/kg) to be released. Such compounds, produced 
by dehydration of solubilized sugars, might serve as nutrient 
sources (Lopez et al. 2004). 


In the current study, the numbers of 16S rRNA gene copies 
increased for RWS1 and TWS1, whereas the numbers remained 
low in the WS control (Fig. 2a). Thus, the bacterial community in 
WS was not stimulated to growth and biodegradation, yet may 
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RWS1 TWS1 


■ Enterobacteriales - Gammaproteobacteria 

■ Pseudomonadales - Gammaproteobacteria 

■ Aeromonadales - Gammaproteobacteria 

■ Flavobacteriales - Bacteroidetes 

■ Sphingobacteriales - Bacteroidetes 

■ Actinomycetales - Actinobacteria 

■ Fungi 

Fig. 4 Venn diagram indicating unique and common bacterial and fungal 
strains in RWS1 and TWS1. Superscripted letter a indicates Raoultellal 
Klebsiella types. Superscripted letter b belongs to Basidiomycota 

have undergone residual growth due to the occurrence of traces 
of nutrients coming from the soil suspension. Overall, the data 
suggested that the lignocellulosic substrate, as compared to the 
WS control, did support bacterial growth (Figs. 1 and 2). With 
respect to fungal communities, we observed a lower capacity to 
respond quickly to the substrates (Fig. 2b). Bacterial abundances 
have been determined earlier using the 16S rRNA gene as a 
proxy, e.g., in lignin-amended soils (DeAngelis et al. 2011a). 
Also, cellulose-degrading bacteria have been quantified using 
qPCR of glycoside hydrolase genes (family GH5 and GH48) 
(Pereyra et al. 2010). Collectively, such methods provide estima¬ 
tions of target gene abundances, yet are subjected to constraints 
as to the sensitivity and specificity related to the primer systems 
used and/or the DNA extraction methods. An assessment of the 
16S rRNA gene copy numbers per cell indicated, on average, 


2.92 and 1.71 copies of 16S rRNA genes per cell for RWS1 and 
TWS1, respectively. It is possible that the presence of a taxon 
such as Bacteroidetes , with an often low 16S rRNA copy 
number (3.5±1.5 copies per genome) compared with 
Gammaproteobacteria (5.8±2.8 copies per genome), strongly 
affected the measured average copy numbers per cell in the 
TWS1 treatment (Vetrovsky and Baldrian 2013). 

PCR-DGGE, on the basis of the 16S rRNA gene and the 
ITS1 region, has been used before on lignocellulose- 
degrading communities (Ho et al. 2012; Ma et al. 2013). 
Recently, Li et al. (2013) reported the use of primers for the 
(3-glucosidase family GH1-3 to evaluate the microbial com¬ 
munity in cattle manure-rice straw composting systems. In 
our study, the initial (inoculum) richness, indicated by the 
band numbers in the bacterial community fingerprints, was 
high (>70), after which divergent patterns emerged for RWS1 
and TWS1, with around 40 % of similarity (based on band 
composition) between these (Fig. 3a). Thus, clearly, the bac¬ 
terial community compositions were driven by substrate treat¬ 
ment. For RWS1, the bacterial communities reached approx¬ 
imate stability after six transfers, whereas for TWS1 only 
minor structural changes were found, related to band numbers 
and positions, after three transfers (Fig. 3a). Similar results 
were reported before (Trifonova et al. 2008), i.e., structural 
stability was found after five transfers in torrefied grass mi¬ 
crobial enrichments. Such dynamics might be due to the 
presence of furanic compounds. A dominant DGGE band 
(B4) from TWS1 was assigned to R. terrigena (Fig. 3a). 
Considering this taxon, which is possibly involved in the 
biodegradation, recent phylogenetic studies demonstrated het¬ 
erogeneity within the genus Klebsiella , suggesting its division 
into the genera Raoultella and Klebsiella (Drancourt et al. 
2001). Thus, band B4 indicated the presence of a Klebsiella- 
like organism, which was dominant and possibly active in 
response to the both wheat substrates. 

With respect to the fungal community structures, the PCR- 
DGGE results showed high community richness in the 



Table 2 Taxonomic affiliation by analysis of partial 18S rRNA gene sequence and enzymatic activity of seven fungal strains recovered in RWS1 and 
TWS1 


Strain code 

Closest relative species—“best” hit partial 18S rRNA 

Phylum 

“Best” hit 
accession 

Identity 

Activity 
in CMC 

Activity in 
xylan 

8wlF a 

Penicillium camemberti strain ATCC 10387 

Ascomycota 

GQ458039 

1,111/1,115(99 %) 

(+) 

(+) 

10t2F b 

Penicillium citrinum 


GU944576 

647/647 (100 %) 

(+) 

(+) 

2wlF a 

Coniochaeta ligniaria (Lecythophora ) strain CBS 267.33T 


AJ496246 

655/660 (99 %) 

(+) 

(+) 

2t2.1F b 

Coniochaeta ligniaria (Lecythophora ) strain CBS 267.33T 


AJ496246 

656/660 (99 %) 

(+) 

(+) 

2t2.2F b 

Plectosphaerella cucumerina strain NRRL 20430 


AF176951.1 

633/643 (99 %) 

(+) 

(+) 

2w2F a 

Plectosphaerella sp. MH727 


FJ430715.1 

649/654 (99 %) 

(+) 

(+) 

2tlF b 

Trichosporon porosum strain CBS 5597 

Basidiomycota 

JN939405 

770/770 (100 %) 

(-) 

(-) 


a Strains recovered from RWS1 
b Strains recovered from TWS1 


4?) Springer 










2798 


Appl Microbiol Biotechnol (2014) 98:2789-2803 



b) 

e 

•g 



■| 




□ Oxidoreductase 
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□ Hemi-cellulase 
■ Total 
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Fig. 5 a Functional screening of 124 bacterial strains recovered from 
RWS1 and TWS1 to detect (hemi)cellulolytic activity (CMC-ase and 
xylanase), growth in mineral medium (MMA) (without-substrate con¬ 
trol), plus glucose (glu), fructose (fru) and furanic compounds (furfural 
and 5-HMF), and oxidoreductase activity in the presence of 
5-HMF by iodine oxidation method ( gray scale color plates). Labels— 
green and red : positive and negative (hemi)cellulolytic activity, 


respectively. Black , gray and white : negative, weak and positive growth, 
respectively. Yellow and blue : positive and negative extracellular oxido¬ 
reductase activity, respectively. Asterisks : negative outside control strains, 
b Number of positive strains with extracellular production of 
hemicellulases, 5-HMF oxidoreductase and ability to grow in 5-HMF in 
each culture (RWS1 and TWS1). Superscripted letter a indicates as 
Raoultella!Klebsiella types 


original soil sample. This richness became markedly reduced communities growing on TWS1 reached stability after the 
throughout the transfers and PC A showed that the second transfer, which was earlier than those from RWS1. 
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The latter showed structural stability after transfer 6 (Fig. 3b). 
The fact that the fungal community structures were very 
dissimilar between RWS1 and TWS1 (-15 % similarity) 
indicated that these were strongly influenced by the substrate 
treatment. The differential stabilization in TWS1 than in 
RWS1 could be explained in the light of the specific conditions 
offered by the substrate, leading to selection of grossly similar 
types fine-tuned to the torrefied material. 

Strains belonging to Serratia , Escherichia , Sphingo- 
bacterium and Arthrobacter were recovered only in the first 
stages of the sequential transfers. Possibly these strains were 
depleted to extinction by the dilution-to-stimulation approach, 
indicating their lack of competitiveness versus other organisms 
in the emerging consortia. Overall, the RWS1 and TWS1 
systems appeared to indeed incite a reshaping of the bacterial 
communities, with reductions in richness and increases in 
prevalence of particular members of the Enterobacteriales , 
Pseudomonadales , Flavobacteriales and Sphingobacteriales 
(Fig. 4 and Table 1). These organisms are likely capable of 
coping with the nutritional and stress constraints of the system, 
similar to what occurs in the mosquito ( Anopheles gambiae ) 
gut microbiome (Wang et al. 2011b). Furanic aldehydes are 
highly reactive molecules, giving rise to the formation of reac¬ 
tive oxygen species (ROS), such as superoxide and H 2 0 2 , 
causing damage to proteins and nucleic acids. Presumably, 
the communities in TWS1 with Enterobacteriales and 
Pseudomonadales may relate to such bacterial stress tolerance. 
With respect to stress (antioxidant) responses, these may be 
primarily coordinated by two regulons (SoxR and OxyR). In 
Enterobacteriales , OxyR is a sensor of H 2 0 2 and it controls the 
activation of major peroxide-degrading enzymes and other 
genes responding to oxidative stress (Kim et al. 2002). Possibly, 
this is the link between the abundance of Enterobacteriales and 
lignin degradation in the RWS1 treatment. Recently, several 
members of the Enterobacteriales were indeed found to pos¬ 
sess genes for plant polymer degradation. They were also found 
in high proportions in insect herbivore microbiomes (Suen et al. 
2010; Aylward et al. 2012). 

We isolated three different ftmgal genera from both systems 
{Coniochaeta , Plectosphaerella and Penicillium ) (Fig. 4 and 
Table 2). The genus Coniochaeta encompasses filamentous 
fungi that inhabit soil and decaying wood. Furanic compound 
utilization (Lopez et al. 2004; Trifonova et al. 2008) and hemi- 
cellulose degradation (Lopez et al. 2007) have been reported for 
such organisms. Recently, Plectosphaerella —which is highly 
related with Acremonium and Sarocladium —has been charac¬ 
terized as a xylanolytic and oleaginous producer yeast 
(Summerbell et al. 2011; Li et al. 2012). Metataxonomic anal¬ 
yses in bagasse feedstock piles showed that the most abundant 
fungal genera were Penicillium , Scytalidiu , Humicola and 
Thermomyces (Rattanachomsri et al. 2011). Eichorst and Kuske 
(2012) performed an elegant experiment, using 13 C maize- 
cellulose, to evaluate the active microbial communities in 


different soils. They identified Trichocladium , Chaetomium, 
Dactylaria and Arthrobotrys as the active fungi. Also, 
Coprinopsis and Pseudallescheria were identified in a micro¬ 
bial consortium growing on wheat straw (Wang et al. 2011a). 
Collectively, such studies point to an important finding, i.e., that 
the microbial source and substrate are highly relevant for the 
ftmgal lignocellulolytic community composition in multispecies 
microbial consortia. Interestingly, we obtained one strain of 
Trichosporon (Basidiomycota) from treatment TWS1 (Fig. 4 
and Table 2). This triglycerol-lipid-producing yeast has been 
shown to be resistant to 5-HMF and furfural (Chen et al. 2009). 
Remarkably, the growth of oleaginous yeasts for lipid produc¬ 
tion using lignocellulose biomass provides a practical option for 
a process with high economic competitiveness (Liu et al. 2013). 

Obtaining strains from a microbial consortium using tradi¬ 
tional plating techniques is restricted to the readily cultivable 
ones and so key components of the consortia might have been 
missed in our work. We performed PCR-DGGE analyses of a 
representative of each BOX-PCR group side-by-side with the 
enriched samples to identify those members of the community 
brought into culture. However, some strains ( Enterobacteriales ) 
produced more than one band in DGGE and the same pattern 
was found between Entembacter , RaoultellalKlebsiella and 
Kluyvera strains, possibly due to high similarity in G+C content 
and identity. However, from the current work, we infer that the 
plating technique used was useful to recover, possibly, the most 
abundant microbial players. 

In the case of the bacteria, 43 and 11 strains, recovered from 
RWS1 and TWS1, respectively, showed (hemi)cellulolytic ac¬ 
tivity (Fig. 5), indicating that the TWS1 system was a rather 
poor selector of such degradative organisms, possibly as a 
result of the presence of furanic compounds, lignin monomers 
and/or oligosaccharides. In methodological terms, Romano 
et al. (2013) reported a correlation between the diameter of 
the degradation halo around colonies and the quantity of re¬ 
ducing sugars produced. However, we used this method as a 
qualitative measurement. 

De Angelis et al. (2011b) recently sequenced the genome 
of Enterobacter lignolyticus SCF1, a novel facultatively an¬ 
aerobic bacterium with the capability to attack lignin. The 
genome showed an arsenal of genes for lignocellulolytic and 
carbohydrate-active enzymes, next to a large number of oli¬ 
gosaccharide transporters. In addition, the genome of 
Flavobacterium johnsoniae has shown an arsenal of enzymes 
for digestion of plant cell wall polysaccharides and hemicel- 
luloses, as identified in silico and experimentally. These re¬ 
sults indicate the ability of these organisms to degrade lignin, 
chitin, starch and cellulose (McBride et al. 2009). Recently, 
Chandra et al. (2011) reported a ligninolytic microbial con¬ 
sortium, composed of Serratia , Klebsiella and Citrobacter 
types, which was capable of decolorizing and detoxifying 
black liquor from the rayon grade pulp manufacturing paper 
industry. Moreover, Ho et al. (2012) showed bio-hydrogen 
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Fig. 6 Diagram for hypothetic bioconversion flux in RWS1 and TWS1. 
Red, green, blue, orange, and white circles : strains belonging to 
Entewbacteriales, Pseudomonadales, Flavobacteriales, Sphingomonadales 
and other orders, respectively. Black circles’, strains belonging to fungi. 
Abbreviations: Acinetobacter (Ac), Aeromonas (Ae), Arthrobacter (Ar), 


Citrobacter (Ci), Coniochaeta (Cn), Entewbacter (En), Escherichia (Es), 
Flavobacterium (FI), Kluyvera (Kl), Penicillium (Pn), Plectosphaerella 
(PI), Pseudomonas (Ps), RaoultellatKlebsiella (Ra), Serratia (Se), 
Sphingobacterium (Sp) and Trichosponon (Tr) 


production using xylan and paper filter as carbon sources by a 
microbial consortium composed of Clostridium , Klebsiella , 
Acinetobacter , Bacillus , Pseudomonas and Ruminococcus 
retrieved from an anaerobic sludge digester. It is important 
to note that Klebsiella is a facultatively anaerobic nonmotile 
bacterium, producing 2,3-butanediol as a major end product of 
the fermentation of glucose (Ji et al. 2011). Additionally, 
Acinetobacter and Pseudomonas strains recovered by mixed 
enrichment cultures with guaiacol-alkali lignin and reed straw 
are able to produce lignin peroxidases (Wang et al. 2013). All 
of these studies highlight the relevance of the microbial 
players coexisting in multispecies microbial consortia. 

Overall, a total of 96 strains retrieved from the two strate¬ 
gies showed growth in 5-HMF (Fig. 5), which exceeded the 
furfural degradation; this suggested that the furfural concen¬ 
tration used was comparatively more toxic than 5-HMF. Sev¬ 
eral species of Cupriavidus , Pseudomonas , Acinetobacter , 
Flavobacterium and Stenotrophomonas are known to de¬ 
grade furanic compounds (Lopez et al. 2004; Wierckx et al. 
2011). Furanic compound degradation may be limited to the 


Gram-negative (aerobic) bacteria, with a few exceptions, such 
as Arthrobacter and Leifsonia spp. We expected a major 
number of bacterial strains (recovered from TWS1) with the 
capacity to grow in 5-HMF. However, the microbial players 
were similar in the two enrichment strategies and possibly 
organisms that are active on lignocellulose have enzymes that 
allow growth on furanic compounds. The metabolism of 5- 
HMF is performed with enzymes encoded by eight so-called 
hmf genes (Wierckx et al. 2011). First, 5-HMF is oxidized to 
the corresponding monocarboxylic acid, which is further ox¬ 
idized to 2,5-furan dicarboxylic acid (FDCA) by the specific 
oxidoreductase HmfH or by nonspecific dehydrogenases 
(Koopman et al. 2010a). This latter transformation entails 
two oxidations and 5-formyl-2-furoic acid may be expected 
as an intermediate product and H 2 0 2 as byproduct (Wierckx 
et al. 2011). Thus, 47 (48.9 %) bacterial strains showed 
extracellular 5-HMF oxidoreductase activity that was revealed 
by an iodide oxidation method (Figs. 5 and 6). According to 
our knowledge, this is the first report of a novel, easy and fast 
method to detect oxidoreductase activity in the presence of 


<£) Springer 






























Appl Microbiol Biotechnol (2014) 98:2789-2803 


2801 


5-HMF. This is important because the 5-HMF oxidoreductase 
might have several functions in the biological oxidation of 
iodide, related with lignin degradation, as well as the produc¬ 
tion of FDCA (new biomass-derived compound that may 
serve as a substitute for terephthalate) (Koopman et al. 
2010b). In the same way, extracellular oxidoreductases can 
detoxify xenobiotics, such as phenolic or anilinic compounds, 
and can efficiently oxidize inorganic compounds such as 
iodide, using oxygen as the final electron acceptor (Karigar 
and Rao 2011). 

To depict the process flux and organisms involved in our 
enrichments, we suggest a hypothetic bioconversion flux 
(Fig. 6). In the RWS1 treatment, the major microorganisms 
involved in the degradation of lignin might be affiliated with 
Pseudomonas , Acinetobacter , Enterobacter and Penicillium . 
Additionally, the (hemi)cellulolytic activity might be related 
with Raoultella/Klebsiella , Kluyvera , Flavobacterium and 
Coniochaeta . Finally, Arthrobacter and Citrobacter could act 
as scavengers of the resulting monomers. In the TWS1 system, 
the (hemi)cellulolytic activity may have been carried out by 
organisms related to Raoultella/Klebsiella and Pseudomonas 
species, whereas the role of Sphingobacterium may be 
scavenger of monomers (Fig. 6). Thus, possibly, mini¬ 
mal consortia functioning in lignocellulose and furanic 
compound transformations could be composed of some 
of the here-isolated organisms which are affiliated with the 
bacterial genera Pseudomonas , Acinetobacter , Raoultella/ 
Klebsiella , Flavobacterium , plus the fungi Coniochaeta and 
Trichosporon . 
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